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Abstract 

Background: Vaccines against HPV16/18 are approved for use in females and males but most countries currently have 
female-only programs. Cultural and geographic factors associated with HPV vaccine uptake might also influence sexual 
partner choice; this might impact post-vaccination outcomes. Our aims were to examine the population-level impact of 
adding males to HPV vaccination programs if factors influencing vaccine uptake also influence partner choice, and 
additionally to quantify how this changes the post-vaccination distribution of disease between subgroups, using incident 
infections as the outcome measure. 

Methods: A dynamic model simulated vaccination of pre-adolescents in two scenarios: 1) vaccine uptake was correlated 
with factors which also affect sexual partner choice ("correlated"); 2) vaccine uptake was unrelated to these factors 
("unrelated"). Coverage and degree of heterogeneity in uptake were informed by observed data from Australia and the USA. 
Population impact was examined via the effect on incident HPV16 infections. The rate ratio for post-vaccination incident 
HPV16 in the lowest compared to the highest coverage subgroup (RR L ) was calculated to quantify between-group 
differences in outcomes. 

Results: The population-level incremental impact of adding males was lower if vaccine uptake was "correlated", however 
the difference in population-level impact was extremely small (<1%) in the Australia and USA scenarios, even under the 
conservative and extreme assumption that subgroups according to coverage did not mix at all sexually. At the subgroup 
level, "correlated" female-only vaccination resulted in RR L = 1.9 (Australia) and 1.5 (USA) in females, and RR L = 1.5 and 1.3 in 
males. "Correlated" both-sex vaccination increased RR L to 4.2 and 2.1 in females and 3.9 and 2.0 in males in the Australia 
and USA scenarios respectively. 

Conclusions: The population-level incremental impact of male vaccination is unlikely to be substantially impacted by 
feasible levels of heterogeneity in uptake. However, these findings emphasize the continuing importance of prioritizing 
high coverage across all groups in HPV vaccination programs in terms of achieving equality of outcomes. 
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Introduction 

Vaccination of pre-adolescent females against human papillo- 
mavirus (HPV) has been recommended or included in publicly 
funded programs in many developed countries. To date, four 
counties have also recommended that pre-adolescent males also be 
vaccinated - Australia, Austria, Canada and the United States [1— 
4], Thus far only Australia has implemented a national publicly 
funded vaccination program which includes boys (from 20 1 3), and 
Austria has announced a publicly funded program including boys 
will commence from 2014 [5]. Many other developed countries 
are still to make decisions around whether or not to fund or 
recommend male vaccination. 

Previous modelling studies have shown that the incremental 
impact of vaccinating males depends on uptake in females, with 
the incremental benefit (and thus cost-effectiveness) decreasing 



with increasing female coverage [6-8], and that increasing 
coverage further in females can be more effective and cost- 
effective than including males in vaccination programs [8-10]. 
The underlying mechanism is that increasing female coverage will 
increase the proportion of heterosexual partnerships where at least 
one partner is vaccinated [8], and so increase indirect protection 
(herd effects). Indirect protection plays an important role in the 
incremental effectiveness of male vaccination compared to female- 
only programs, because as female coverage increases, so too does 
indirect protection for males and unvaccinated females; therefore 
the extent to which vaccinating males can offer additional 
protection becomes progressively smaller. However, if factors 
influencing vaccine uptake also influence sexual partner choice, 
this would potentially affect the extent of indirect protection, and 
as a result alter the incremental impact of male vaccination. To 
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date, the large majority of modelling studies have assumed that 
choice of sexual partner and sexual behavior are not correlated 
with vaccination uptake - that is, that any factors which may be 
associated with vaccine uptake are not also associated with aspects 
of sexual behavior. However, this assumption is unlikely to hold in 
practice since sociodemographic, cultural and geographic factors 
are all likely to have a substantial bearing on both sexual partner 
choice and vaccination uptake in some settings. Such associations, 
may be, but are not necessarily, mediated by levels of sexual 
activity. For example, geographic factors are likely to impact on 
partner choice across different behavioral groups and also are 
likely to play a role in vaccine coverage, particularly where 
vaccination delivery is coordinated at a regional or state level. 

Two recent systematic reviews of factors associated with vaccine 
uptake in teenage girls, based primarily on studies conducted in 
the USA, identified race, area of residence, cost/insurance status, 
physician's recommendation, completion of other recommended 
vaccinations, and parental concerns about vaccine safety as 
important factors [11,12]. These factors have also been identified 
in more recent studies [13-18]. Factors associated with uptake 
may vary depending on the method of delivery (school- versus 
primary care/ clinic-based), funding model (publicly-funded versus 
insurance- or private-payer) and between countries. For example, 
a recent study from Canada found that socioeconomic status was 
associated with uptake where the vaccine delivery was clinic-based, 
but not when it was school-based [19]. This is also supported by 
data from two settings with school-based publicly-funded pro- 
grams which found little effect of socioeconomic status on uptake 
[20,21]; in contrast, in the United States public funding is limited 
and insurance status appears to affect uptake [11-13]. The extent 
to which ethnicity is associated with uptake has varied between 
studies [11,12], and the underlying reason behind any identified 
associations with ethnicity may differ between countries. For 
example, a recent review for the United Kingdom found some 
evidence of variation in acceptability and uptake by ethnicity, but 
that this may in practice have been a proxy for religion [22]. Data 
from the USA suggest that non-religious factors underpin the 
association with ethnicity, as parental intention to vaccinate their 
child did not vary by ethnicity, whereas uptake did [18]. Since 
male HPV vaccination has not yet been widely adopted 
internationally, studies looking at factors associated with uptake 
in young males have been very limited. However, associations are 
likely to be similar to those observed for females [23,24], and 
factors such as race, cost/insurance status, physician's recommen- 
dation and completion of other recommended vaccinations have 
also emerged in two recent studies of uptake in young males in the 
USA [14,25]. Parents whose daughters had already received the 
HPV vaccine were more willing for their sons to receive it [26,27] 
and gender of child is not generally associated with acceptability to 
parents [28,29]. Therefore, it is plausible that factors influencing 
uptake in females are likely to also influence uptake in males. 

The simplifying assumption that factors influencing vaccine 
uptake are not also associated with aspects of sexual behavior may 
have impacted the accuracy of prior estimates of both the 
effectiveness of female vaccination and the incremental benefit of 
male vaccination. No prior study has performed a comprehensive 
evaluation of the impact of correlation between vaccine uptake 
and sexual behaviour (whether driven by cultural, geographic or 
other factors) on estimated vaccination outcomes. Two prior 
Canadian studies have assessed the direct impact of differential 
vaccination uptake in sexual behavior subgroups (i.e. subgroups 
defined according to rates of partner change and age of onset of 
sexual activity) for female [30] and both-sex [31] vaccination. 
Malagon el al [31] found that disparities in uptake would 



potentially reduce population-level effectiveness if uptake was 
lower in more sexually active groups, but that this effect 
decreased at higher levels of vaccine coverage. However, no 
studies to date have considered the impact of factors influencing 
partner choice other than, as above, the direct impact of 
differential coverage according to absolute level of sexual activity 
(number of sexual partnerships), and nor has any study assessed 
the extent to which including males affects the differences in 
outcomes between groups in a program with differences in female 
uptake. Additionally, quantifying between-group differences has 
rarely been done in modelling studies; these have generally 
focused on population-level effectiveness and cost-effectiveness, 
because these are critical factors in policy decisions. In this 
regard, however, a highly influential report to the WHO has 
recommended that health equity impact be routinely assessed for 
proposed policies [32]. 

The aims of the current study were, therefore, to perform a 
generalized assessment of the impact of a correlation between 
factors influencing vaccine uptake and choice of sexual partner; 
where such a correlation could be underpinned by a number of 
plausible mechanisms driving changes in both outcomes, including 
geographic, socioeconomic and cultural elements. We examined 
the impact at the population level of adding males to HPV 
vaccination programs if factors influencing vaccine uptake also 
influence partner choice. Additionally, we aimed to quantify how 
these specific kinds of differential uptake may change the 
distribution of disease between subgroups after vaccination, using 
incident infections as the outcome measure. Inequalities between 
groups here are described, and represent between-group differ- 
ences in either vaccination uptake or infection outcomes after 
vaccination. However it should be noted that inequality is not by 
itself directly interpretable as "inequity", whichis often seen as 
inequality which is avoidable and unfair, [33]. 

Methods 

Model description 

A dynamic model of HPV 16 was used to simulate sexual 
behavior and HPV vaccination, transmission and natural history 
in a population. This model has been previously described in detail 
[7,34,35]. Briefly, the model simulates the potential transmission 
of HPV during heterosexual partnerships in a population which is 
closed and stratified by sex, age and level of sexual activity. Model 
parameters were obtained via literature reviews and fitting to 
observed data [7,34]. 

Simulations were performed under various conditions, to reflect 
either scenarios where vaccine uptake within the population was 
correlated with factors which also affect choice of sexual partners 
("correlated"), or where uptake was unrelated to any of these 
factors ("unrelated"). In the "correlated" scenario, subgroups were 
simulated in separate closed models which did not interact 
sexually, to represent mixing which is assortative with respect to 
factors correlated with vaccine uptake. Different levels of vaccine 
uptake were applied in different subgroups, but all other model 
parameters were consistent between subgroup populations (for 
example sexual behavior, including the proportion of the 
population in different sexual behavior strata, and all aspects of 
vaccine efficacy and HPV natural history). In the corresponding 
"unrelated" scenario the equivalent population-level vaccine 
uptake, which was derived by weighting the subgroup coverage 
rates by the respective sizes of the subgroups, was applied across 
the population. 
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Vaccination coverage scenarios 

This analysis does not explicitly model any specific factors 
associated with vaccine uptake; rather, this was intended to be an 
exploratory example from which broad conclusions could be 
drawn. This was done because it appeared that factors associated 
with vaccine uptake (encompassing both the aspects of initiation 
and completion of the three-dose series) vary in different settings. 

However, the scenarios chosen were informed by data from 
actual HPV vaccination experiences. Three population coverage 
scenarios were considered: "Australia" (72% overall); "USA" 
(32% overall); and a hypothetical exploratory scenario with 
intermediate coverage (50% overall) where the "correlated" 
variant explored "extreme inequality". The "Australia" coverage 
scenario was based on observed vaccine uptake in Australia for 
girls aged 15 years in 2009 (who were offered vaccination in 2007 
and 2008). Three subgroups were formed, with their size and 
vaccine coverage based on three-dose uptake within the different 
states, grouped based on broad coverage levels [36,37]. Uptake in 
the subgroups ranged from 64.5% to 76.3%. The "USA" 
population coverage scenario was based on observed vaccine 
uptake for females in the USA who were aged 13 to 17 years in 
2010. Four subgroups were formed, with their size and vaccine 
coverage likewise based on three-dose uptake within the different 
states grouped based on broad coverage levels [38,39]. Uptake in 
the subgroups ranged from 23.0% to 43.6%. In each coverage 
scenario, the equivalent overall population coverage used in the 
"unrelated" variant of the scenario was the weighted average of 
three-dose coverage across all subgroups. Detailed information on 
coverage assumptions in each case is presented in Figure 1 and 
Table SI in File SI. 

In each case, we compared results from a female-only program 
to those from a both-sex program. We also examined the extent to 
which including males in a program altered the distribution of 
disease between groups in the population after vaccination. All 
scenarios assumed that vaccination was delivered at age 1 2 years. 
In the baseline analysis, we assumed equal coverage in males and 
females, but assessed lower coverage in males relative to females in 
a sensitivity analysis. 




igh 



Lower coverage Higher coverage Moderate coverage, 
scenario ("USA") scenario Extreme inequality 

("Australia") 



Figure 1. Subgroup size and vaccine uptake in modelled 
coverage scenarios. Values next to bar represent coverage in that 
subgroup; bar height represents subgroup size. 
doi:1 0.1 371 /journal.pone.01 01 048.g001 



Assessing outcomes and quantifying between-group 
differences in outcomes 

The age-standardized rate (ASR) of incident HPV 16 infections 
at post-vaccination equilibrium (i.e. several decades after the 
implementation of vaccination) was used as the outcome measure 
and a measure of post-vaccination disease burden. Outcomes from 
the model were examined at both the population level and at the 
subgroup level. Outcomes at the population level were calculated 
for the "correlated" uptake scenario by weighting the outcomes 
from die individual subgroups according to their relative size in 
the population. Comparisons were then made with outcomes from 
the corresponding "unrelated" uptake scenario. 

For comparison of outcomes in different subgroups, two 
approaches were used. In the first, the cumulative proportion of 
the population represented by each subgroup was plotted against 
the cumulative proportion of the burden of disease experienced by 
that subgroup (pseudo Lorenz curve). As previously, incident 
HPV 16 infections at post-vaccination equilibrium were used as the 
measure of post-vaccination disease burden; these infection rates 
were assumed to be equivalent in all subgroups (but not in all 
sexual behavioral, sex and age strata within subgroups) prior to 
vaccination. In a population where outcomes are equal across the 
subgroups, the plotted pseudo Lorenz curve would be a diagonal 
line; the further away this plot is from the diagonal line of equality, 
the greater the degree of inequality in outcomes. The distance 
between the ideal diagonal line and the situation which arises in 
each scenario can be quantified by calculating a pseudo Gini 
coefficient, a measure of dispersion in outcomes between different 
groups, which is equal to twice the area between the plotted 
pseudo Lorenz curve and the ideal diagonal [40]. A coefficient of 
zero represents perfect equality, whereas the closer the value is to 
the theoretical maximum, the more unequal outcomes are in 
subgroups. Since outcomes are considered at the level of a 
subgroup (rather than the individual level), the theoretical 
maximum of the coefficient will depend on the subgroup sizes. 
Here the theoretical maximum was calculated to represent the 
most unequal situation possible, which was assumed to occur if all 
disease occurred in the subgroup with the lowest vaccine uptake. 
The Lorenz curve and Gini coefficient are well recognized 
measures of dispersion, and the pseudo Gini coefficient has been 
extensively used in measuring health inequalities [40-42]. 

The second approach to comparing outcomes between 
subgroups calculated the rate ratio (RRtJ as the ratio of the 
post-vaccination equilibrium age-standardized rate of incident 
HPV 16 infections in the subgroup with the lowest vaccine 
coverage relative to that in the subgroup with the highest 
coverage. We additionally used a rate ratio because this is a 
commonly understood concept in health; however it is a less 
desirable measure of inequality, as it does not take into account the 
whole population, nor the sizes of the subgroups with differing 
outcomes [40], 

Sensitivity analysis 

Parameters which were found to have the most influence on 
outcomes in prior work [7], including vaccine coverage and 
aspects of vaccine efficacy, were varied during sensitivity analysis. 
Natural history assumptions were not varied, as previous work 
showed that using alternative model parameter sets which had 
been fitted to HPV prevalence had little impact on predicted 
vaccination outcomes [7]. Details of the assumptions used for 
sensitivity analysis are available in the accompanying Supplemen- 
tary Material (see Sensitivity Analyses in File SI). Additional 
sensitivity analyses varied the extent of heterogeneity in vaccine 
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Figure 2. Impact of heterogeneity in vaccine uptake on population level outcomes. (A) Female-only program (50% overall coverage, 
extreme inequality). (B) Both sex program (50% overall coverage, extreme inequality). "Correlated" uptake refers to a situation where vaccine uptake 
within the population is correlated with factors which also affect choice of sexual partners. "Unrelated" uptake refers to a situation where vaccine 
uptake is unrelated to any of these factors. Vaccination was assumed to commence in 2007. 
doi:1 0.1 371 /journal.pone.01 01 048.g002 



uptake in the "Australia" and "USA" coverage scenarios, in either 
both sexes or males only (see Exploratory Analyses in File SI). 

Results 

Population level outcomes 

In a female-only program with 50% overall coverage, the 
predicted long-term relative reductions in female and male HPV 
infections are 56% and 49% for "correlated" uptake under 
assumptions of 'extreme inequality', compared to 62% and 41% 
for "unrelated" uptake, respectively (Figure 2a). If males were also 
vaccinated at the same coverage as females, the predicted 
reductions in female and male infections are 61% and 60% for 
"correlated" uptake with extreme inequality, compared to 79% 
and 78% for "unrelated" uptake, respectively (Figure 2b). Similar 
effects were seen in the "Australia" (higher) and "USA" (lower) 
coverage scenarios, where the heterogeneity was less pronounced, 
but in these scenarios the effects were extremely small (< 1 % 
difference in population impact) (Table 1). Generally "correlated" 
uptake resulted in a lower vaccination effectiveness at the 
population level, except that female-only programs with "corre- 
lated" uptake resulted in better outcomes in males (but not in 



females), than female-only programs with equivalent "unrelated" 
uptake. The extent to which this occurred varied depending on the 
degree of heterogeneity in uptake between subgroups, but can be 
explained by considering the implications in the simple "extreme 
inequality" scenario. In this example, the assumption that 
coverage occurred in equal-sized subgroups, one with higher 
(90%) and one with lower (10%) coverage resulted in a net loss of 
benefit for females, because the extent of indirect (herd) protection 
was lower in both of the subgroups than in the scenario with 
moderate (50%) coverage in females which was uniform across the 
population. Females in the lower coverage subgroup did not 
experience a large herd effect because coverage was very low; and 
within the higher coverage subgroup, there was less room for 
improvement in outcomes due to herd effects, since coverage, and 
impact, was already very high for females. In contrast, males who 
partnered within the subgroup of females with very high uptake 
experienced strong herd effects. As the benefits to males from a 
female-only program do not scale linearly with coverage, and 
appear to accelerate as coverage in females increases [43,44], the 
gains in this subgroup of males outweighed the loss for males in the 
lower coverage subgroup. 
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Figure 3. Impact of heterogeneity of vaccine uptake on subgroup outcomes. (A) Higher population coverage ("Australia"; 72.4% overall). (B) 
Lower population coverage ("USA"; 32.1% overall). (C) 50% overall coverage, extreme inequality. 
doi:1 0.1 371 /journal.pone.01 01 048.g003 




Figure 4. Distribution of disease outcomes (incident HPV16 infections) across subgroups (pseudo Lorenz curve). (A) Higher 
population coverage ("Australia"; 72.4% overall). (B) Lower population coverage ("USA"; 32.1% overall). (C) 50% overall coverage, extreme inequality. 
Comparison of the proportion of disease borne by each subgroup with the group's size. The diagonal line represents a situation where there are no 
inequalities in outcomes between subgroups; the further away a plot of outcomes is from this equality line, the more unequal outcomes are in that 
scenario. The pseudo Gini coefficient represents twice the area between the pseudo Lorenz curve and the equality line. 
doi:1 0.1 371 /journal.pone.01 01 048.g004 
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Figure 5. Impact of varying model assumptions on inequality 
of outcomes. (A) Higher population coverage ("Australia"; 72.4% 
overall). (B) Lower population coverage ("USA"; 32.1% overall). A higher 
value of the pseudo Gini coefficient represents more unequal 
outcomes. (F) denotes the value for the pseudo Gini coefficient relating 
to outcomes in females; (M) denotes the value for the pseudo Gini 
coefficient relating to outcomes in males. SA= sensitivity analysis. * 
Switched heterogeneity: Higher heterogeneity used for Australia 
scenario ((equivalent to heterogeneity in main USA scenario); Lower 
heterogeneity used for USA scenario (equivalent to heterogeneity in 
main Australia scenario). 
doi:1 0.1 371 /journal.pone.01 01 048.g005 



Subgroup level outcomes 

Compared to female-only vaccination, including males in the 
vaccination program improved outcomes in all subgroups, as the 
absolute risk of HPV 1 6 infection reduced in all subgroups and in 
both sexes (Figure 3). 

As expected, female-only vaccination with "correlated" vaccine 
uptake resulted in differential outcomes between subgroups, and 
this was more pronounced for outcomes in females than in males 
(Figure 4). However, in the higher and lower coverage scenarios 
examined, the differences in outcomes between subgroups were 
comparatively small (Table 1). Including males in the program 
further increased the differences between subgroups after vacci- 
nation, and the level of inequality became similar for males and 
females (Figure 4). The absolute value of the pseudo Gini 
coefficient increased for males and females when males were 
included in programs, indicating outcomes had become more 
unequal between the subgroups, both in females and males 
(Table 1). The extent of this varied, however, and was more 
pronounced in the higher "Australia" coverage scenario than the 
lower "USA" coverage scenario, as in the higher coverage 
scenario the post-vaccination incidence of disease in the subgroup 
with highest coverage was extremely small (Figure 5). In the lower 
coverage scenario, the pseudo Gini coefficient increased from 
0.0771 to 0.1258 for females and from 0.0439 to 0.1200 for males 
when males were also vaccinated (theoretical maximum for four 



subgroups modelled: 0.8378; Table 1). In the higher coverage 
scenario, the pseudo Gini coefficient increased from 0.0936 to 
0.2205 for females and from 0.0578 to 0.2086 for males when 
males were vaccinated (theoretical maximum for three subgroups 
modelled: 0.8766). In the scenario with moderate coverage but 
extreme inequality, the pseudo Gini coefficient increased from 
0.4696 to 0.4997 for females and 0.4002 to 0.4996 for males when 
males were vaccinated (theoretical maximum for two subgroups 
modelled: 0.5). 

Using an alternative measure of inequality, the rate ratio, 
produced similar results; that is, that "correlated" female 
vaccination led to differential subgroup outcomes in females and 
also to a lesser extent in males, and that including males increased 
the between-group differences in outcomes. Female-only vaccina- 
tion with "correlated" vaccine uptake resulted in RR L =1.9 
("Australia" scenario), 1 .5 ("USA" scenario) or 31.9 (extreme 
inequality scenario) in females, and RRl= 1.5,1.3 and 9.0 in these 
scenarios in males. Both-sex vaccination with "correlated" uptake 
increased RR L to 4.2, 2.1 and 3,317 in females and 3.9, 2.0 and 
2,467 in males in programs with higher, lower, or extreme 
inequality in coverage respectively (Table 1). 

Sensitivity analysis 

The predicted reductions in female and male infections, and the 
incremental benefit achieved by adding males, remained very 
similar at the population level for "correlated" and "unrelated" 
uptake in the "Australia" and "USA" scenarios during sensitivity 
analysis. The difference in the predicted percentage reduction in 
infections at the population-level from a program with "unrelated" 
versus "correlated" uptake was less than one percentage point in 
almost all cases. The exception was a version of the "Australia" 
scenario where heterogeneity in vaccine uptake was greater 
(equivalent to that observed in USA data), although the difference 
was still quite small (less than five percentage points). Sensitivity 
analysis demonstrated that adding males to female-only programs 
consistently improved the absolute outcomes in all subgroups, 
regardless of overall coverage level, or duration of vaccine 
protection. The incremental benefit of adding males was greatest 
when vaccine protection was long-lasting in males, but not in 
females. However, compared to the female-only program, adding 
males also consistendy increased between-group differences in 
outcomes, although the extent of this effect varied (Figure 5). 
Scenarios where coverage was lower in males, or duration of 
vaccine protection was short in males, had less impact on the 
degree of inequality. However the reason for this was because in 
practice these programs were less effective over the long term, and 
so the relatively greater benefits experienced by the subgroups with 
higher uptake were short-lived and did not increase between- 
group differences over the long term. Adding males with uniform 
uptake also had less impact on the degree of inequality, however 
this was still increased compared to female-only vaccination 
(Figure 5; see also File SI). 

The magnitude of between-group differences in outcomes was 
influenced both by the extent of inequality in vaccine uptake, and 
also by the overall population coverage level. For programs with 
equivalent heterogeneity in vaccine uptake, between-group 
differences in outcomes were greater in a program with higher 
coverage than one with lower coverage (Figure 5; also see 
Exploratory Analyses in File SI). 

Discussion 

Some factors which have high potential to influence HPV 
vaccine uptake include geographical location within a country 
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(since programs are often organized at a regional or state level), 
cultural identity (since programs may be more or less successful at 
reaching different groups) and sociodemographic factors (since, 
particularly when vaccination delivery is not school based, 
sociodemographic factors may affect the accessibility of or decision 
to be vaccinated) [11,12,19,21,36,38,45]. These factors are also 
likely to influence the social networks in which people congregate 
and thus the pool from which potential sexual partners are chosen. 
The impact of a partial correlation between sexual partner choice 
and vaccine uptake on vaccination effectiveness has not previously 
been assessed. We found that inequalities in vaccine uptake which 
are correlated with factors that also affect choice of sexual partner 
do not adversely impact the population level effectiveness of 
vaccination programs, unless the inequalities in uptake are 
pronounced. Reassuringly, this finding suggests that the estimated 
cost-effectiveness of both female and male HPV vaccination, based 
on evaluations which to date have assumed homogenous uptake, is 
unlikely, except in extreme situations, to be substantially affected 
by this kind of inequality. Although "correlated" uptake generally 
resulted in a worse outcome at the population level, the difference 
was extremely small in most scenarios we explored. The 
comparatively small differences at the population level can be 
explained by the comparatively small degree of heterogeneity in 
uptake in the observed data from Australia and the USA, and thus 
the comparatively small differences which arose in the probability 
that at least one partner in a heterosexual relationship was 
vaccinated. Herd effects are related to the chance that at least one 
partner in a relationship is vaccinated [8], and the degree of 
heterogeneity modeled in our "Australia" and "USA" coverage 
scenarios did not substantially alter the overall probability that at 
least one partner would be vaccinated in a heterosexual 
partnership compared to the "unrelated" scenario. Interestingly, 
female-only programs with "correlated" uptake resulted in better 
population-level outcomes in males (but worse population-level 
outcomes in females), than female-only programs with equivalent 
"unrelated" uptake. In effect, "correlated" uptake in a female-only 
program appears to shift some of the vaccination benefits from 
females to males (relative to those which would be achieved by a 
"unrelated" program), specifically, shifting benefits to the male 
partners of females in the subgroups with higher coverage. 

Our study confirmed that vaccinating males improved out- 
comes in all subgroups, even if uptake was correlated with factors 
that also affect choice of sexual partner. However, we also found 
that in this situation including males in vaccination programs 
tended to further concentrate disease in particular subgroups and 
increase between-group differences in outcomes, even if uptake is 
uniform in males. In these situations, decision-makers may face a 
trade-off between reducing inequalities versus reducing the levels 
of disease overall. It would likely then be important to take into 
account the highly setting-specific factors of who the subgroups 
are, their relative advantage and disadvantage overall, their other 
risk factors, and the potential alternative uses of resources, and 
whether the inequality is avoidable or not. That is, decision makers 
would need to take into account whether or not the inequalities 
represent inequity in a particular setting, given that inequity is 
often and additionally defined as inequality which is avoidable and 
unfair [33]. For example if infections were likely to become more 
concentrated in subgroups that were already disadvantaged, or in 
females less likely to attend cervical screening and thus at higher 
risk of cervical cancer, including males may be seen as less 
favorable than alternative interventions (including, for example, 
interventions to achieve higher uptake in the disadvantaged 
subgroups). Conversely, reducing HPV infections further in a 
particular subgroup could be beneficial if individuals in that group 



are otherwise disadvantaged, or less likely to be screened. We 
could not and have not examined inequity per se in this 
generalized assessment, because this would require consideration 
of setting-specific issues. Setting-specific analyses could incorporate 
measures which have been proposed to quantify inequity, rather 
than inequality, such as the concentration index, which could not 
be used here as they require subgroups to be ordered according to 
a gradient of socioeconomic status or advantage [40] . 

We have used the post-vaccination rate of new HPV infections 
as a proxy for disease risk in different subgroups, since it provides 
an outcome which is comparable for males and females. The 
current analysis does not simulate long term outcomes of 
infections, such as cancer. This is because the risk of developing 
an HPV-related cancer depends on a range of further factors; in 
particular cervical cancer risk would be strongly influenced by 
cervical screening behavior. The specific impact on cervical cancer 
would also depend on whether variations in uptake are correlated 
with screening uptake, and the direction of the relationship. For 
example in the United States, HPV vaccine uptake is higher in 
areas where screening coverage is higher [46], whereas in New 
Zealand some ethnic groups with lower screening coverage have 
higher HPV vaccine uptake [45,47], We have also only modelled 
one HPV type, HPV16, in this analysis. HPV16 is associated with 
the highest risk of oncogenesis in humans and this approach is 
consistent with other exploratory model-based analyses [10], and 
done on the basis that analyses for other types would give 
qualitatively similar results. 

A previous study by Malagon et al [31] assessed the population- 
level impact of disparities in vaccine uptake if vaccine uptake 
differed according to level of sexual activity (i.e. number of 
partnerships and age at initiation), whereas we performed a 
complementary assessment of the effects of differences in vaccine 
uptake according to any factor that influences the choice of specific 
partner(s). Our results were similar in finding that pronounced 
heterogeneities in vaccine uptake may affect the population-level 
impact of HPV vaccination programs. Our finding that including 
males tended to increase inequality relative to a female-only 
program differed from the findings of the previous study; however 
this particular finding from Malagon et al [31] was based on a 
scenario with uniform uptake across subgroups where inequalities 
were driven by varying herd effects in the subgroups resulting from 
their different levels of sexual behaviour (since even uniform 
uptake in females increased inequality). While there were 
differences in the precise behavioural subgroups modelled and 
factors associated with vaccine uptake between the two studies, an 
exploratory analysis suggested that another reason for the different 
findings was that different measures were used to quantify 
inequality. One measured the absolute level and distribution of 
disease between groups, whereas the other measured the relative 
impact of vaccination. Information about the exploratory analysis 
and a more detailed discussion of these issues are available in 
Supplementary Material (see Exploratory Analyses in File SI). 

An argument in favor of male vaccination is that it would 
improve equality between males and females. In order to fully 
explore equality in relation to sex, it would be necessary to 
consider additional disease endpoints, such as HPV-related 
cancer, genital warts and recurrent respiratory papillomatosis, to 
fully capture the differences in the burden of HPV-related disease 
by sex. While we have not focused on the differences in outcomes 
in males versus those in females here, several points are worth 
considering. Firstly, modelling studies have consistently shown that 
targeting one sex for HPV vaccination is more effective at 
reducing the overall burden of disease in the population (as well as 
being more cost-effective) [8-10]. Secondly, the burden of HPV- 



PLOS ONE | www.plosone.org 



8 



August 2014 | Volume 9 | Issue 8 | e101048 



Impact of Inequalities in HPV Vaccine Coverage 



related disease is generally much greater in females than in males, 
and this is why vaccination of females only (rather than males only) 
is more effective and cost-effective [10,48]. Thirdly, the obvious 
inequality arising from a female-only vaccination program is not in 
males generally (since their burden of HPV-related disease was 
lower than that of females prior to vaccination, and is likely to 
decrease as a result of female vaccination due to herd effects 
[7,8,10,49]); rather, it will specifically be in men who have sex with 
men (MSM). MSM have a higher burden of HPV-related disease 
than other males, particularly anal disease [50,51], and are less 
likely to benefit from herd effects from a female-only program 
[52]. Vaccination of MSM has been shown to be cost-effective in 
at least one setting, even in early adulthood when targeted 
vaccination may be more feasible [53,54]. This is an important 
topic for further research in settings with female-only vaccination, 
although local factors such as the disease burden and HIV 
prevalence will be influential [53]. 

In our model, the subgroups we considered were assumed not 
to interact sexually. This did not mean that vaccinated 
individuals only partnered with other vaccinated individuals (or 
unvaccinated only with unvaccinated), or that sexual behavior 
was homogenous within the subgroups; instead it implied that 
subgroups with varying vaccine uptake (potentially representing 
geographic, cultural or socioeconomic subgroups) did not interact 
sexually. In reality, subgroups are not closed, however our 
modelling assumption allowed us to explore the opposite extreme 
assumption to that previous models of male vaccination have 
generally made, which is that uptake is uniform, or equivalently 
that factors associated with vaccine uptake are unrelated to 
factors associated with choice of partner. Where inequalities in 
uptake were not very great, population-level effectiveness was 
very similar to the "unrelated" variant, even under our extreme 
and conservative simplifying assumption that subgroups did not 
intermix, and so this result should apply to the more realistic 
situation where subgroups are not completely closed. However, 
given that the assumption of closed subgroups was an extreme 
assumption for our main population-level question, between- 
group differences will be more muted in the more realistic 
situation where subgroups are not closed. 

In conclusion, this analysis suggests that the incremental 
impact (and therefore the cost-effectiveness) of female-only or 
male vaccination will not be substantially impacted by hetero- 

References 

1. Centers for Disease Control and Prevention (201 1) Recommendations on the use 
of quadrivalent human papillomavirus vaeeine in males-Advisory Committee on 
Immunization Practices (ACIP), 2011. MMWR MorbMortalWklyRep 60: 
1705-1708. 

2. National Advisory Committee on Immunization (2012) Update On Human 
Papillomavirus (HPV) Vaccines. Can CommunDis Rep 38. 

3. Pharmaceutical Benefits Advisory Committee (2011) November 2011 PBAC 
Outcomes - Positive Recommendations. 

4. European Centre for Disease Prevention and Control (2012) Introduction of 
HPV vaccines in EU countries - an update. Stockholm, Sweden. 

5. Bundesministcrium fur Gesundheit (Austrian Federal Ministry for Health) (2013) 
Ausbau Kinder-Impfprogramm: Die Aufhahme der Impfung gegen Humane 
Papilloma- Viren (HPV) [Expanding children's vaccination program: the 
inclusion of vaccination against human papilloma virus (HPV)] . 

6. Canfcll K, Chcsson H, Kulasingam SL, Bcrkhof J, Diaz M, ct al. (2012) 
Modeling Preventative Strategies against Human Papillomavirus-Related 
Disease in Developed Countries. Vaccine 30 Suppl 5: F157-F167. 

7. Smith MA, Lew JB, Walker RJ, Brotherton JM, Nickson C, ct al. (201 1) The 
predicted impact of HPV vaccination on male infections and male HPV-related 
cancers in Australia. Vaccine 29: 91 12-9122. 

8. Brisson M, van de Velde N, Franco EL, Drolet M, Boily MC (201 1) Incremental 
impact of adding boys to current human papillomavirus vaccination programs: 
role of herd immunity. J Infect Dis 204: 372-376. 



geneity in population uptake unless there is pronounced 
inequality in vaccine uptake that is associated with factors 
influencing partner choice. Outcomes in all subgroups can be 
improved by vaccinating males as well as females, but including 
males could increase the degree of inequality in outcomes 
between population subgroups, whether the variations in uptake 
between subgroups carry through to males or not. Therefore, 
reducing inequalities in vaccine uptake remains an important 
goal, and its importance is not diminished by extending female- 
only vaccination programs to include males, since including 
males may exacerbate rather than mitigate between-group 
differences in outcomes. In general terms these findings 
emphasize the continuing importance of achieving high 
coverage in all geographic, cultural and socioeconomic groups 
in HPV vaccination programs. 

Supporting Information 

File SI Supporting text and tables. Table SI. Subgroup 
size and vaccine uptake in coverage scenarios modelled. Table 
S2. Summary of main results in exploratory analysis, by sex, 
coverage scenario and program type. Additional Analyses 
Performed - Sensitivity Analyses and Exploratory Anal- 
yses. 
(DOCX) 

Acknowledgments 

Wc gratefully acknowledge the contributions of other members of the 
Cancer Modelling Group at UNSW and external collaborators to the 
development of the general model structures used in this evaluation. 

Both authors were previously employed in the Cancer Research 
Division, Cancer Council NSW, Sydney, Australia. 

We also thank the investigators of the Australian Study of Health and 
Relationships and the Australian Social Science Data Archive for provision 
of data on sexual behavior in Australia, which were used to inform aspects 
of the dynamic model used in this paper, and who bear no responsibility for 
the further analysis and interpretation of these data. 

Author Contributions 

Conceived and designed the experiments: MS KC. Performed the 
experiments: MS. Analyzed the data: MS. Wrote the paper: MS KC. 

9. Chcsson HW, Ekwucmc DU, Saraiya M, Dunne EF, Markowitz LE (201 1) The 
cost-effectiveness of male HPV vaccination in the United States. Vaccine 29: 
8443-8450. 

10. BogaardsJA, Krctzschmar M, Xiridou M, Mcijer CJ, Bcrkhof J, ct al. (2011) 
Sex-specific immunization for sexually transmitted infections such as human 
papillomavirus: insights from mathematical models. PLoS Med 8: el 00 1147. 

1 1 . Kcsscls SJ, Marshall HS, Watson M, Braunack-Mayer AJ, Reuzel R, ct al. 
(2012) Factors associated with HPV vaccine uptake in teenage girls: A systematic 
review. Vaccine. 

12. Fisher H, Trotter CL, Audrey S, MacDonald-Wallis K, Hickman M (2013) 
Inequalities in the uptake of human papillomavirus vaccination: a systematic 
review and meta-analysis. Int J Epidemiol 42: 896—908. 

13. Dcmpsey AF, Schaffcr SE, Cohn LM (2012) Follow-up analysis of adolescents 
partially vaccinated against human papillomavirus. J Adolesc Health 50: 421- 
423. 

14. Gilkey MB, Moss JL, McRee A-L, Brewer NT (2012) Do correlates of HPV 
vaccine initiation differ between adolescent boys and girls? Vaccine 30: 5928- 
5934. 

15. Hamlish T, Clarke L, Alexander KA (2012) Barriers to HPV immunization for 
African American adolescent females. Vaccine 30: 6472—6476. 

16. Keenan K, Hipwell A, Stepp S (2012) Race and sexual behavior predict uptake 
of the human papillomavirus vaccine. Health Psychology 31: 31—34. 

17. Liddon NC, Lcichliter JS, Markowitz LE (2012) Human papillomavirus vaccine 
and sexual behavior among adolescent and young women. Am J Prev Med 42: 
44-52. 



PLOS ONE | www.plosone.org 



9 



August 2014 | Volume 9 | Issue 8 | e101048 



Impact of Inequalities in HPV Vaccine Coverage 



18. Sadigh G, Dempsey AF, Ruffin M, Resnicow K, Carlos RC (2012) National 
patterns in human papillomavirus vaccination: an analysis of the National 
Survey of Family Growth. Hum Vaecin Immunothcr 8: 234—242. 

19. Musto R, Siever J, Johnston J, SeidelJ, Rose M, et al. (2013) Social equity in 
Human Papillomavirus vaccination: a natural experiment in Calgary Canada. 
BMC Public Health 13: 640. 

20. Barbaro B, Brotherton JM, Gertig DM (2012) Human papillomavirus 
vaccination and cervical cancer screening by socioeconomic status, Victoria. 
McdJ Aust 196: 445. 

21. Kumar VM, Whynes DK (2011) Explaining variation in the uptake of HPV 
vaccination in England. BMC Public Health 11: 172. 

22. Marlow LA (2011) HPV vaccination among ethnic minorities in the UK: 
knowledge, acceptability and attitudes. Br J Cancer 105: 486^92. 

23. Liddon N, Hood J, Wynn BA, Markowitz LE (2010) Acceptability of human 
papillomavirus vaccine for males: a review of the literature. J Adolesc Health 46: 
113-123. 

24. Ogilvie GS, Remple VP, Marra F, McNeil SA, Naus M, et al. (2008) Intention of 
parents to have male children vaccinated with the human papillomavirus 
vaccine. Sex Transm Infect 84: 318-323. 

25. Reiter PL, Gilkey MB, Brewer NT (2013) HPV vaccination among adolescent 
males: Results from the National Immunization Survey-Teen. Vaccine 31: 
2816-2821. 

26. Reiter PL, McRcc AL, Gottlieb SL, Brewer NT (2010) HPV vaccine for 
adolescent males: acceptability to parents post-vaccine licensure. Vaccine 28: 
6292-6297. 

27. Reiter PL, McRcc AL, Kadis JA, Brewer NT (2011) HPV vaccine and 
adolescent males. Vaccine 29: 5595-5602. 

28. Zhang SK, Pan XF, Wang SM, Yang CX, Gao XH, et al. (2013) Perceptions 
and acceptability of HPV vaccination among parents of young adolescents: A 
multiccnter national survey in China. Vaccine 31: 3244-3249. 

29. Lenselink CH, Gerrits MM, Melchcrs WJ, Massuger LF, van Hamont D, et al. 
(2008) Parental acceptance of Human Papillomavirus vaccines. Eur J Obstet 
Gynecol Reprod Biol 137: 103-107. 

30. Shafcr LA, Jeffrey I, Elias B, Shearer B, Canfell K, et al. (2013) Quantifying the 
impact of dissimilar HPV vaccination uptake among Manitoban school girls by 
ethnicity using a transmission dynamic model. Vaccine 31: 4848-4855. 

3 1 . Malagon T, Joumier V, Boily M-C, Van de Velde N, Drolct M, et al. (20 1 3) The 
impact of differential uptake of HPV vaccine by sexual risks on health 
inequalities: A model-based analysis. Vaccine 31: 1740-1747. 

32. Commission on Social Determinants of Health (2008) Closing the gap in a 
generation: health equity through action on the social determinants of health. 
Final Report of the Commission on Social Determinants of Health. Geneva. 

33. Whitehead M (1992) The concepts and principles of equity and health. 
Int J Health Serv 22: 429-445. 

34. Smith MA, Canfell K, Brotherton JM, Lew JB, Barnabas RV (2008) The 
predicted impact of vaccination on human papillomavirus infections in 
Australia. Int J Cancer 123: 1854-1863. 

35. Walker R, Nickson C, LewJB, Smith M, Canfell K (2012) A revision of sexual 
mixing matrices in models of sexually transmitted infection. StatMed 31: 3419- 
3432. 

36. Department of Health and Ageing (2011) Human Papillomavirus (HPV): 
Immunise Australia program. Department of Health & Ageing. 



37. Australian Bureau of Statistics (2008) Population by Age and Sex, Australian 
States and Territories, Jun 2002 tojun 2007. 

38. Centers for Disease Control and Prevention (2011) National, State and Local 
Area Vaccination Coverage among Adolescents Aged 13 through 17 Years - 
United States 2010. Morbidity- and Mortality Weekly Report 60: 1 1 17-1 123. 

39. U. S. Census Bureau (201 1) Annual Estimates of the Resident Population for the 
United States, Regions, States, and Puerto Rico: April 1, 2010 to July 1, 2011 
(NST-EST20 11-01). 

40. Wagstaff A, Paci P, van Doorslacr E (1991) On the measurement of inequalities 
in health. Soc Sci Med 33: 545-557. 

41. Preston SH, Haines MR, Pamuk E (1981) Effects of Industrialization and 
Urbanization of Mortality in Developed Countries. IUSSP 19th International 
population Conference. Manila: IUSSP. pp. 233-254. 

42. Leclerc A, Lert F, Fabien C (1990) Differential mortality: some comparisons 
between England and Wales, Finland and France, based on inequality measures. 
Int J Epidemiol 19: 1001-1010. 

43. Garnett GP, Kim JJ, French K, Goldie SJ (2006) Chapter 21: Modelling the 
impact of HPV vaccines on cervical cancer and screening programmes. Vaccine 
24 Suppl 3: S3/178-S173/186. 

44. Garnett GP (2005) Role of herd immunity in determining the effect of vaccines 
against sexually transmitted disease. J InfectDis 191 Suppl 1: S97— 106. 

45. Smith L, Grigg M, Davies L, Reilly J, Laukau S (2012) HPV Immunisation 
Programme Implementation Evaluation. Wellington, NZ. 

46. Jemal A, Simard EP, Dorell C, Noone AM, Markowitz LE, et al. (2013) Annual 
Report to the Nation on the Status of Cancer, 1975-2009, Featuring the Burden 
and Trends in Human Papillomavirus (HPV)-Associated Cancers and HPV 
Vaccination Coverage Levels. J Natl Cancer Inst 105: 175-201. 

47. Smith M, Walker R, Canfell K (2013) National Cervical Screening Programme 
Annual Report (2010-2011). Wellington, NZ. 

48. Forman D, de Martel C, Lacey CJ, Socrjomataram I, Lortet-Tieulent J, et al. 
(2012) Global burden of human papillomavirus and related diseases. Vaccine 30 
Suppl 5: F12-23. 

49. Ali H, Donovan B, Wand H, Read TR, Regan DG, et al. (2013) Genital warts in 
young Australians five years into national human papillomavirus vaccination 
programme: national surveillance data. BMJ 346: f2032. 

50. Palefsky JM (2010) Human Papillomavirus-Rclatcd Disease in Men: Not Just a 
Women's Issue. Journal of Adolescent Health 46: S12— S19. 

51. DalingJR, Madeleine MM, Johnson LG, Schwartz SM, Shera KA, et al. (2004) 
Human papillomavirus, smoking, and sexual practices in the etiology of anal 
cancer. Cancer 101: 270-280. 

52. Donovan B, Franklin N, Guy R, Grulich AE, Regan DG, et al. (2011) 
Quadrivalent human papillomavirus vaccination and trends in genital warts in 
Australia: analysis of national sentinel surveillance data. Lancet Infect Dis 11: 
39-44. 

53. Kim JJ (2010) Targeted human papillomavirus vaccination of men who have sex 
with men in the USA: a cost-effectiveness modelling analysis. Lancet Infect Dis 
10: 845-852. 

54. Simathcrai D, Bradshaw CS, Fairley CK, Bush M, Heley S, et al. (2009) What 
men who have sex with men think about the human papillomavirus vaccine. Sex 
Transm Infect 85: 148-149. 



PLOS ONE | www.plosone.org 



10 



August 2014 | Volume 9 | Issue 8 | e101048 



